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ABSTRACT 

We present the measured Sunyaev-Zel'dovich (SZ) flux from 474 optically-selected MaxBCG clusters that 
fall within the Atacama Cosmology Telescope (ACT) Equatorial survey region. The ACT Equatorial region 
used in this analysis covers 510 square degrees and overlaps Stripe 82 of the Sloan Digital Sky Survey. We 
also present the measured SZ flux stacked on 52 X-ray-selected MCXC clusters that fall within the ACT 
Equatorial region and an ACT Southern survey region covering 455 square degrees. We find that the measured 
SZ flux from the X-ray-selected clusters is consistent with expectations. However, we find that the measured 
SZ flux from the optically-selected clusters is both significantly lower than expectations and lower than the 
recovered SZ flux measured by the Planck satellite. Since we find a lower recovered SZ signal than Planck, 
we investigate the possibility that there is a significant offset between the optically-selected brightest cluster 
galaxies (BCGs) and the SZ centers, to which ACT is more sensitive due to its finer resolution. Such offsets 
can arise due to either an intrinsic physical separation between the BCG and the center of the gas concentration 
or from misidentification of the cluster BCG. We find that the entire discrepancy for both ACT and Planck can 
be explained by assuming that the BCGs are offset from the SZ maxima with a uniform random distribution 
between and L5 Mpc. Such large offsets between gas peaks and BCGs for optically-selected cluster samples 
seem unlikely given that we find the physical separation between BCGs and X-ray peaks for an X-ray-selected 
subsample of MaxBCG clusters to have a much narrower distribution that peaks within 0.2 Mpc. It is possible 
that other effects are lowering the ACT and Planck signals by the same amount, with offsets between BCGs 
and SZ peaks explaining the remaining difference between ACT and Planck measurements. Several effects that 
can lower the SZ signal equally for both ACT and Planck, but not explain the difference in measured signals, 
include a larger percentage of false detections in the MaxBCG sample, a lower normalization of the mass- 
richness relation, radio or infrared galaxy contamination of the SZ flux, and a low intrinsic SZ signal. In the 
latter two cases, the effects would need to be preferentially more significant in the optically-selected MaxBCG 
sample than in the MCXC X-ray sample. 

Subject headings: cosmic microwave background- galaxies: clusters: general - galaxies: clusters: intracluster 
medium 
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1. INTRODUCTION 

Galaxy cluster properties may follow simple scaling laws 
reflecting their self-similarity. This possibility has given 
credence to their use as cosmological probes. Cosmolog- 
ical parameters have been obtained by X-ra y, optical, and 
most recently SZ surveys of clusters (e.g., 'Vikhlin in et al.l 
■^09; Mantz et al. 2010; Rozo et alJl201 0: Sehgal et al.l l2011l: 
Benson et al.ll201 II) . At the same time the scaling laws that 



feed into these parameter constraints continue to undergo 
scrutiny. 

Recent millimeter-wavelength data have opened a new 
window whereby these scaling relations can be robustly 
checked against SZ flux measurements. The SZ cluster sig- 
nal has been predicted to have a low-scatter correlation with 
cluster mass (e.g., iMotl et al.1 120051: iNagaii 2006). If true, 
this would make SZ-detected clusters an excellent tracer of 
structure growth in the Universe (e.g ., Wang & Steinhardt 
T998HHaiman et alJlIOOlt [Holder et alBoOlMCarlstrom et al.i 



20021) . As steps towards understanding the SZ-mass rela- 
tion, several studies have shown that the SZ clust er sig- 



nal correlates well with X -ray sign als (e.g., Bonamen te et alJ 
12008, 2011; Andersson et al.ll20I ll) . dynamically determined 
masses CSifon et al. 201^7 and weak-lensing determined 
masses (e.g..lMarrone et al1 l201 ll) . 

In particular, the Planck satellite recently reported a 
good agreement between the measured and expected SZ- 
mass relation for a sarnple of X-ray-selected clusters 
(iPlanck Collaboration et al.l 1201 lallbl) . However, a similar 
comparison for optically-selected clusters yielded an am- 
plitude of SZ flux lower than expected by about a fac- 
tor of t wo, with an even larger discrepanc y for lower-mass 
clusters (Planck Collaboration et al."201 IcV An analysis by 
[Dra per et al. (201 1 ) using data from the WMAP satellite found 
a simi lar result, however with larger uncertainty. Ha nd et alJ 
(l20Tll) . using data from t he Atacama Cosmology Telescope 
(ACT) jSwetz et alj|201 ll) and stacking luminous red galax- 
ies, also suggested a low SZ flux for optically-selected halos. 
Among the possible explanations could be that either the SZ 
signal is not a robust tracer of galaxy clusters and groups or 
that optical selection techniques are somehow biased. 

Here, we investigate this discrepancy by stacking optically- 
selected clusters in millimeter- wavelength data fro m ACT that 
overlaps Stripe 82 of the Sloan Digital Sky Survey (lYork et alj 
|M)0). We also measure the SZ flux for X-ray-selected clus- 
ters as a consistency check. Understanding these scaling rela- 
tions will have important implications for cluster astrophysics 
as well as for their use in cosmological studies. 

This paper is organized as follows. Section |2]discusses the 
data sets used in this analysis. Section[3]describes the method 
used to measure cluster SZ flux. Results are presented in Sec- 
tions|4]through|6]and discussed in Section]?] 

2. DATA SETS 

Below we describe the catalogs of optically-selected and X- 
ray-selected clusters and the millimeter-wavelength data used 
to measure the cluster SZ fluxes. We note that throughout this 
work Msooe refers to the mass within Rsooc, which is the ra- 
dius within which the average density equals 500 times the 
critical density of the Universe at the cluster redshift. Sim- 
ilarly, M200H1 gives the mass within R2Q0m, the radius within 
which the average density equals 200 times the mean mat- 
ter density of the Universe at the cluster redshift. A fidu- 
cial cosmology of ftm = 0.27, ft^ = 0.73, and h = 0.71 is 



also adopted (iKomatsu et alJl20Tll) . where H(z) = HoE(z) = 
(h X 100 km s-i Mpc-i)£(z) and E(z) = [flM +z? + Qa]^^^. 

2.1. The MaxBCG Optical Cluster Catalog 

The MaxBCG Optical Cluster Catalog consists of 13,823 
clusters selected from Data Rele ase 5 (DR5) of the Sloan Dig- 
ital Sky Survey (iKoester et al. l 2007a b). The clusters were 
selected from a 7500 deg^ area of sky using the observation 
that cluster galaxies tend to be the brightest galaxies at a given 
redshift, share a similar red color, and are spatially clustered. 
The catalog consists of clusters that fall in the redshift range 
of 0.1 < z < 0.3 and have a richness measure, A^zoom, within 
10 < A^zoom < 190. The richness is defined as the number of 
red-sequence galaxies with L > 0.4L* (in the / band) within 
a projected radius of Rimm- The catalog provides the BCG 
position (RA and DEC), photometric redshift, richness, BCG 
luminosity, and total luminosity of each cluster. Applying the 
cluster detection method to mock catalogs suggests that the 
catalog should be 90% pure and 85% complete. 

Mass estimat es of the clus ters in the MaxBCG sample 
were d erived by iSheldon et al.l (2009) and Mandelbaum et al.] 
( 2008a ) usin g weak gravitation al lensing. Johnston et al] 
(l2007al) and iRozo et al] (120091) used those mass determi- 
nations to c onstru ct richness-mass (Af200m -^sooe) relations. 
iRozo e t aT (2009), in particular, used the masses derived 
from iMandelbaum et al. (2008a) due in part to the au- 
thors' c areful treat ment of p hotometric r e dshi ft uncertainties 
( IMandelbaum et alj|2008bl) . IRozo et alj (12009 ) also stacked 
the MaxBCG clu ster catalog on X-r ay maps from the ROSAT 
AU-Sky Survey (Voges et al."1999^ and used the Lx-M re- 
lation from Vikhlinin et al. (2009) as a prior to inform their 
richness-m ass relation. Thus the richness-mass relation of 
IRozo et al] (i200 9) is expected to be consistent w ith an Lx-M 
relation from X-ray clusters. IRozo et aD (12010b applied this 
richness-mass relation to the MaxBCG sample of optically- 
selected clusters and found a cosmological constraint on erg of 
as(n,„/0.25)°-^^ = 0.832 ± 0.033 assuming a flat ACDM cos- 
mology. Throughout this work we define t he Nypom -M^nor re - 
lation as given by Eqs. 4, A20, and A21 of IRozo et all (120091) . 



2.2. The MCXC X-ray Cluster Catalog 

The Meta-Catalog of X-ray detected Clusters of galaxies 
(MCXC) is presented in iPiffaretti et al.l (.2011). The MCXC 
cluster catalog is based on pubUcly available data from a num- 
ber of different X-ray catalogs including the ROSAT AU-Sky 
Survey and comprises 1743 clusters. The catalog provides 
the position (RA and DEC), redshift. X-ray 0.1-2.4 keV 
band luminosity (Lsooe), mass (Msooc), and radius (Rsooc) of 
each system. The redshift distribution of this catalog goes 
from about 0.05 to 1. The Lx-M relation derived from the 
MCXC clusters in ^ffaretti et al"] (120111) is co nsistent with 
fliat of Pratt et al." (2009^) and IVikhhnin et al] (|2009). The 
l^khlinin et al. (2009) Lx-M relation was used to derive a 
cosmology constraint on erg from a sample of X-ray clusters 
that is a subsample of the MCXC catalog. From this analy- 
sis di ey found o-8(^»,/0.25) ° '*^ = 0.813 ±0.013 (stat) ±0.024 
(sys) (I Vikhhnin et al.l 120091) . Other authors have found sim- 
ilar cons traints on as frorn ROSAT and othe r X-ray cluster 
samples dHenrv et al.l2009HMantz et alJiMol) . 

2.3. Millimeter-wave Data from the Planck Satellite 

Given that the optically-selected MaxBCG cluster catalog 
and the X-ray-selected MCXC cluster catalog yield consistent 
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Lx -M relations and constraints on trs, one would expect both 
cluster samples to yield consistent Jsoot- - A^socv relations PI 
However, in a set of papers presented by the Planck collabo- 
ration (Planck Collaboration et al. 201 Ib c) it was found that 
the F500C ~ ■'WsoOf relation for X-ray-selected clusters from the 
MCXC sample agrees with expectations, whereas the normal- 
ization of the Fgooc - A^sooc relation for the optically-selected 
MaxBCG sample was lower than expectations by about a fac- 
tor of two. For both cases, e xpectations were ba sed on X-ray 
derived cluster profiles from Arnau d et al.l (120101) . For the lat- 
ter case,Jhe expectation also folded in the A^200m -^500e rela- 
tion of lRozo et aij (120091) . The Planck data used in the above 
analysis consists of six HFl c hannel millimeter-wave tem- 
peratu re maps as described in iPlanck HFI Core Team et aP 
( I2OI Ih . This data set comprises the first ten months of the 
survey and covers the full sky. 

2.4. Millimeter-wave Data from the Atacama Cosmology 
Telescope 

The Atacama Cosmology Telescope is a six-meter tele- 
scope operating at an altitude of 5200 meters in the Atacama 
Desert of Chile. The telescope site allows ACT to observe in 
both the northern and southern hemispheres. In this work, we 
use millimeter-wave maps covering two regions of sky: one 
spanning 510 deg^ over the celestial equator and one span- 
ning 455 deg^ in the southern hemisphere. The Equatorial 
region consists of a 4.5°-wide strip centered at a declina- 
tion of 0° and running from 20''20"' through O'' to 03''50'". 
The Southern region consists of a 7° wide strip centered on 
-53° and extending from 00''12"' to 7''10"'. Both sky regions 
were observed over the 2008, 2009, and 2010 observing sea- 
sons at 148 and 218 GHz. The Equatorial region overlaps the 
Sloan Digital Sky Survey (SDSS) Stripe 82 and thus overlaps 
492 clusters in the MaxBCG catalog. The Equatorial plus 
Southern regions combined overlap 74 clusters in the MCXC 
catalog For a more detailed description of the ACT in- 
strument, observations, an d data reduction see Fowler et aD 
_itz et al.' ("201 W 'Marriage e t al.l f201 li) : ,Das et ail 
Hajianet al. (201 1); Dunner et al.1 (120121) . 
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3. MEASUREMENTS OF SZ FLUX 

3.1. Multi-frequency Matched Filter 

We use a multi-frequency matched filter to extract 
the thermal SZ sig nal from clusters as described in 
iHaehnelt & Tegmarklrfl996h and lMelin et aPdlOOeh . The fil- 
ter in Fourier space is given by 

*(k) = a2 [P(k)r'.T(k), 
where T(k) has the components 



(1) 



T,(k) = T'(k);V B,(k). 



(2) 



Here is the frequency dependence of the thermal SZ sig- 
nal for frequency v, t'{V) is the profile of the cluster in 
Fourier space, and Bi/(k) is the profile of the instrument beam 
in Fourier space. P(k) is the power spectrum of the noise, 
both astrophysical and instrumental. The astrophysical noise 

J'sooi- is the SZ flux within Rsoor. dividing out the frequency dependence 
of the SZ signal. 

Only 6 clusters are in common between the 492 cluster MaxBCG sample 
and the 74 cluster MCXC sample. 



sources for cluster detection include the primary lensed mi- 
crowave background, radio galaxies, infrared galaxies. Galac- 
tic emission, and the SZ background from unresolved clus- 
ters, groups, and the intergalactic medium. Since the power 
from the cluster thermal SZ signal is subdominant to these 
astrophysical sources (as evide nced by Luekeret al. (2 01(1); 
iHall et all (l207ol) ; iFowler et al.l (12010); Das etiD (1201 ll) ). we 
approximate the power spectrum of the total noise as the 
power spectrum of the data itself. Here 



1 



(fk [T(k)]'-[P(k)]->-[T(k)] 



(3) 



.(2^)2 

is the normalization of the filter that ensures an unbiased esti- 
mate of the cluster signal. 

3.2. SZ Model Template 

We use for the filter's spatial template the empirical 
universal pressure profile of Arnaud et al. (2010) derived 
from X-ray observati ons of the REXCESS cluster sample 
dBohringer et al.ll2007l) . The three-dimensional pressure pro- 
file is given by 



1 



jc'^(l+jc")(^-T)/" 



(4) 



where x= r/r,, r^ = Rsooc/csoo, C500 = 1.156, a = 1.0620, /3 = 
5.4807, and 7 = 0.3292. The normalization of this profile is 
arbitrary for the purposes of the matched filter We essentially 
measure this normalization for each cluster when we apply 
this filter to our maps. The SZ signal is given by the projected 
gas pressure, so we describe the filter template by integrating 
the three-dimensional profile above along the line of sight. 
Thus 



p2O(0)= / 2P^''{^Jf + e^DA(zf)dl, (5) 
Jo 

where l^ax = 5^500c and Da{z) is the angular diameter distance. 
The filter is truncated at 5R5ooc/Da(z) = 5 0^00^, which contain s 
over 95% of the signal, as was don e in iMelin et aP (1201 ll) : 
IPlanck Collaboration etaTI (1201 Ibid) . 

3.3. Application of Filter 

Before filtering the maps, we establish uniform noise prop- 
erties by creating an effective weight map that has pixel-wise 
effective weights given by Weff = {■^+ ■^)~' where w, is the 
pixel weight for the ith frequency. Here weight is defined as 
the number of observations per pixel normalized by the obser- 
vations per pixel in the deepest part of the map. We multiply 
the 148 and 218 GHz ACT data maps pixel-wise by the square 
root of the effective weight map. After we apply the filter to 
create a filtered map, we divide the filtered map pixel-wise by 
the square root of the effective weight map. 

We appl y the matched filter above follow ing the procedure 
given in Pl anck Collaboration et al.l (1201 Id) . For each cluster 
in the MaxBCG catalog that falls in the ACT coverage region, 
we create a unique matched filter using the A^200m - sooc re- 
lation of Rozo et al. (2009) to determine each cluster's M500C 
and subsequently Rsqoc from the A^20()m value given in the cat- 
alog. We also derive Da(z) for each cluster from its photo- 
metric redshift. When the filter is applied to the map, the 
pixel coincident with the location of the cluster center in the 



4 



Sehgal et al. 



G 10-2 



o 

Si 
CO 



10-3 



o 
a 
S 

o 
o 
in 



10- 



10- 



10- 



—\ 1 T" 



• ACT Sim 
■ Model 



10 



100 



N200 



Fig. 1. — Recovered yjoo values from 447 simulated clusters embedded in 
ACT 148 and 218 GHz maps (black circles). These clusters were simulated 
to mat ch the properties of the clusters in the MaxBCG catalog I Koester et al. 
I2007af) that overlap the ACT equatorial region. The simulated clusters were 
placed at random locations within the ACT maps. The input F500 values are 
shown as blue squares. 



filtered map should have a value equal to the normalization 
of the two-dimensional SZ template given by Eq. |5] To sim- 
plify extraction of the desired quantity, we normalize P^'^iO) 
itself to equal unity when integrated over 9 in two dimensions 
from zero to 56*5001.. Thus the pixel value recovered at the 

C\'l 

cluster center position after applying the filter is F^g . Here 



■■ y^j^^^^ , where Y^j^^^^ is the integrated projected 
SZ signal within a cylinder of radius 5/?5ooc- We use a geomet- 
ric factor of y^fot = (0.986 / 1 .8 ^WIr,^^. given in Appendix A 
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Fig. 2. — The black solid Une shows the distribution of offsets between 
BCG and X-ray gas peak from 208 clusters found in both the MaxBCG and 
MCXC cluster catalogs. The red dashed line shows the distribution for the 
subsample of rich clusters (N2oom > 35). The blue dotted line shows the 
same for the subsample of poor clusters (N2Q0m < 35). There are 100 and 108 
clusters in each subsample respectively. 
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Fig. 3. — Black circles show recovered $500 values from 446 embedded 
simulated clusters where the assumed cluster centers used for Y^qq recovery 
are offset from the input centers with a random distribution given by the black 
solid line in Figure [5] "Mcent" stands for miscentered. The input values are 
shown as blue squares. 



of iMelin et al] (1201 Ih to convert from Y"^[ 



to ¥''''' , the in- 



tegrated SZ flux within a sphere of radius /?5ooc- Through- 
out this work we plot ?5oo = F5oo£"2/^(z)(ZJa(z)/500 Mpc)2, 
where isoo = Yg^^^, as in lPlanck Collaboration et alJ (1201 Id) . 

4. SIMULATED ACT SZ SIGNALS 

4. 1 . Embedding Simulated Clusters in ACT Maps 

In order to test the analysis pipeline which applies the SZ 
extraction procedure discussed above, we use simulated clus- 
ters embedded within the ACT data maps at random locations. 
Using the information in the MaxBCG catalog for the 492 
clusters that fall within the ACT Equatorial region, we create 
a unique SZ profile for each cluster using Eq. |5]above and the 
cluster Rsooc and z given in the catalog. We then add each 
simulated SZ cluster to the 148 and 218 GHz ACT Equatorial 
maps, placing it at a random location and scaling the thermal 
SZ signal to give it the appropriate frequency dependence in 
each map. The simulated SZ signal is also convolved with the 
appropriate ACT beam prior to embedding it within each ACT 
map. We then exclude any simulated clusters from further 
analysis that happen to be within 5' of a point source detected 
at greater than in either the 148 or 2 1 8 GHz maps. We also 
exclude any clusters that are within 10' from the edge of the 
map. As an additional cut, we exclude all clusters that are in 
noisy parts of the map where the local value of the effective 
weight map is less than 15% of the maximum value. These 
cuts leave 447 clusters. Similarly below, when we extract the 
SZ signal from the real cluster positions using the MaxBCG 
and MCXC catalogs, we apply the same cuts discussed above. 

In Figure [T] we show the results of this SZ extraction pro- 
cedure. Here the simulated clusters are binned by the opti- 
cal richness given in the MaxBCG catalog. The blue squares 
show the input model, and the black circles show the recov- 
ered signal of the simulated clusters embedded at random lo- 
cations within the ACT maps. The error bars are given by 
a/^/N where is the variance of fsoo in each richness bin, 
and is the number of clusters in each richness bin. The vari- 
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Fig. 4. — Measure d Fjoo values for 52 MCXC X-ray-selected clusters 
jPiffaretti et al.ll20Tl1) that fall within the ACT equatorial and southern sur- 
vey regions (black circles). Also shown are expected F500 values based on 
meas ured cluster X-ray properties (blue squares). A cluster profile model 
fromlAmaud et al. (2010) was assumed for determining both measured and 
expected 7500 values. 



ance dominates the uncertainty in each bin, as demonstrated 
in Figure 4 of Planc k Collaboration et al. (2011c ) There is 
good agreement between input and recovered signals. 

4.2. Effect of Cluster Miscentering 

The exercise above shows we should expect excellent re- 
covery of the SZ flux for clusters given in the MaxBCG cata- 
log provided the A^200m -A^sooe and Msooc-isooc relations are 
correct, and the cluster properties (position, redshift, A^200m) 
listed in the catalog a re accurate. However, o ne source of un- 
certainty identified in iJohnston et al.l ( | 2007b) (section 4.3) is 
the positional accuracy of the cluster center. In the MaxBCG 
catalog, the cluster position is given by the location o f the 
Brightest Cluster Galaxy (BCG). ljohnston et all (l2007bl) sug- 
gested two reasons why the BCG found by the MaxBCG clus- 
ter identification algorithm may be offset from the true clus- 
ter center One is that the true BCG may be intrinsically 
offset from the dark matter center or center of the gas con- 
centration presumably due to unrelaxed behavior (e.g., merg- 
ers). Another reason may be that the BCG could b e misiden- 
tified by the cluster finder IJohnston et al.l (12007 b) explore 
this latter effect with mock optical cluster catalogs and find 
that a richness-dependent fraction of the clusters have accu- 
rately identified BCGs while the rest are miscentered due to 
BCG misidentification following a Rayleigh distribution with 
a scale parameter a equal to cruise = 0.42/z"'Mpc. The dis- 
tribution of the intrinsic BCG offset from the gas center is 
unknown|3 

We explore the potential offset between BCG and gas cen- 
ter by studying the subset of clusters in common to both the 
MaxBCG and MCXC catalogs. From the fuU catalogs (which 
have 13,823 and 1743 clusters respectively), we identify 208 

^" Since no scatter in the A'200m ~ ^500r relation has been included here, 
these error bars reflect only the SZ flux recovery error. 

Note that this offset does not arise from pointing uncertainties in optical, 
X-ray, or millimeter-wave instruments. It is due to either BCG misidentifica- 
tion or cluster astrophysics. 
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Fig. 5. — Measured Y^qq values for 474 MaxBCG optically-selected clus- 
ters that fall within the ACT equatorial survey region (black circles). Ex- 
pected 7500 values are shown as blue squares. Both measured and ex- 
pected values assume the A^oOm ~ ^500c relation from Rozo et al.. (2009|) 
and the Arnaud et al. (2010) cluster profile. Red triangles are the measured 
values from the Planck satellite f or a sample of 13,104 MaxBCG clusters 
IPlanck CoUaboration et alJ2011cll . 



clusters that are in common in both catalogs. Here we define 
a cluster as matched in both catalogs when the identified clus- 
ter redshifts are within Az < 0.015 and the projected physi- 
cal separation of the identified cluster centers is less than 1 .5 
Mpc. Using these 208 clusters, we plot the fraction of clus- 
ters as a function of separation between BCG and X-ray peak 
in Figure |2] The solid black line shows the offset distribu- 
tion for the 208 clusters. Half of the clusters have offsets less 
than 0. 1 Mpc, while the other half have a roughly flat offset 
distribution between 0.1 and 1.5 Mpc PI 

The 208 clusters were then divided into rich clusters (with 
A^20()/>i > 35) and poor clusters (with A^200/ii < 35). A richness 
cut of A^2()0m = 35 divides the 208 clusters into roughly even 
subsamples of ^ 100 clusters each. The dashed red line in 
Figure|2]shows the offset distribution for the rich clusters, and 
the dotted blue line shows the distribution for the poor clus- 
ters. Similar distributions are found for both subsamples, with 
the poor clusters showing slightly more offset. 

Using the simulations described above, we explore the ef- 
fect of cluster miscentering on SZ flux recovery. We use the 
same simulated clusters embedded at random positions in the 
148 and 218 GHz ACT maps as before. However, when re- 
covering the SZ fluxes, we use positions for the clusters that 
differ from the true cluster gas centers with a random distri- 
bution given by the black solid line in Figure|2] We also allow 
for 10% of the clusters to be false d etections to ma tch the pu- 
rity of the MaxBCG catalog ( Koester et alj|2007bi) . Figure [3] 
shows the result of the SZ recovery process. The black circles 
show the recovered isooc values for simulated clusters embed- 
ded in ACT data. To give a sense of the relation between Mpc 

We note that the choice of 1.5 Mpc is somewhat arbitrary. When we 
allow matches within 1 Mpc, we find 189 clusters. We also find that if we 
allow matches within 3 Mpc, Figure |2] plateaus instead of dropping to zero 
at large separation. The plateau is due to poor clusters, as the distribution 
of rich ones does tend to zero. The size of a cluster (R20Q111) is usually less 
than 2 Mpc, so this suggests that some of these poor cluster matches may be 
spurious. 
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Fig. 6. — The black open circles show the measured Yjoo values from simu- 
lations with an offset between identified cluster center and gas concentration 
center given by a uniform random distribution between and 1 .5 Mpc. The 
solid black circles, blue squares, and red triangles are the same as in Figure 
\5\ The values from simulations have been shifted in the x-axis for claiity. 

and arcminutes, for the redshift range of the MaxBCG clus- 
ter sample (z e (0.1,0.3)), 0.5 Mpc corresponds to about 2' to 
4.5'. Given the ACT beam size (1.4' at 148 GHz and 1' at 218 
GHz), we expect a decrease in the recovered Fjooc signal due 
to the amount of miscentering shown in Figure |3] How much 
the recovered SZ signal decreases depends on the noise prop- 
erties of the ACT data, which differ significantly from pure 
white noise. We discuss this further in Section]?] 

5. MEASURED ACT SZ SIGNALS 

5.1. Stacking X-ray Selected Clusters 

Within the 510 deg^ ACT Equatorial region and the 455 
deg^ ACT Southern region are located 74 cluste rs found in 
the MCXC catalog of X-ray-selected clusters (Piff aretti et alj 
|2011). After making the cuts discussed in Section 143] to ex- 
clude clusters near bright point sources, near map edges, or in 
very noisy parts of the map, 52 MCXC clusters remain. Using 
the Rsooc, Msmc, and redshift of each cluster and the projected 
SZ profile given in Eq.15] we calculate the expected mean ysooc 
values in each Msmc bin shown as the blue squares in Figure 
|4] The black circles in Figure |4] show the mean of the re- 
covered Y500C values using the multi-frequency matched filter 
given in Eq. [T]and the projected SZ profile created uniquely 
for each cluster The error bars are the error on the mean 
given by a j\fN, where cr^ is the variance and is the num- 
ber of clusters in each bin. Figure |4]shows overall agreement 
between expected and recovered SZ signals for the X-ray- 
selected clusters. The reduced chi-squared is 0.76 using 7 de- 
gr ees of freedom. This is consistent wi th the agreement found 
by iPlanck Collaboration et aP (1201 Ibh for a larger sample of 
X-ray-selected clusters. 

5.2. Stacking Optically Selected Clusters 

In the ACT Equatorial region there are 492 MaxBCG clus- 
ters. This reduces to 474 clusters once the above men- 
tioned cuts are made|3 Using the A^200hi - Msooc rela- 

Note that the ACT region is centered on Stripe 82, and extends beyond 
it. So fewer real clusters are cut than simulated clusters because the real 



tion of | Rozo et all (l2009t) and the M^qoc - ^sooc relation of 
lArnaud et al.l ( l2010h . we find the expected and recovered Y^ooc 
values of the MaxBCG clusters using the method described in 
Section[3] Figure|5]shows the expected values as blue squares 
and the recovered values as black circles. The recovered sig- 
nal is significantly lower than the expected signal as well as 
the signal reco vered by Planckfiom 13, 104 MaxB CG clusters 
(red triangles) (iPlanck Collaboration et al.ll201 Icb . 

From Figure |3] we see that some amount of offset between 
BCG and gas peak can result in a lower measured signal than 
expected. However, this figure also shows that the offset dis- 
tribution given by Figure does not result in a low enough 
measured signal to explain the measurement shown in Fig- 
ure |5] To investigate the amount of offset necessary to match 
the ACT measured isoo values shown in Figure |5] we redo 
the analysis shown in Figure |3] However, this time we use 
an offset distribution that is uniformly random between and 
1.5 Mpc. We also again allow for 10% false detections. The 
results are shown as the open black circles in Figure |6] The 
solid black circles, blue squares, and red triangles in Figure|6] 
are the same as in Figure |5] We see that this amount of off- 
set between BCG and gas peak roughly matches the measured 
values. An extensive scan of offset distributions is beyond the 
scope of this work, but it may be that a more complex or re- 
fined distribution could give a better fit to the measurements. 

5.3. Stacking Optical Clusters Using BCG Dominant 
Subsample and New Richness Measure 

Using a subsample of clusters with "dominant BCGs" we 
examine whether the expected SZ signal is closer to the mea- 
sured values. Such a subsample may more closely correspond 
to an X-ray -selected subsample, and with such a subsample 
IPlanck Collab oration et al. (201 Ic) found better agreement 
between model and meas urement. We follow the definition of 
"BCG dominant" used by lPlanck Collaboration et al.1 (1201 Id) 
which is defined relative to the quantity LBCo/C^'tot-^BCG)- 
Here Ltot and Lbcg are the R-band luminosities of the cluster 
and cluster BCG respectively. For a "BCG dominant" cluster, 
this ratio is larger than the average ratio for a given richness 
bin. From the sample of 474 MaxBCG clusters used above, 
126 are "BCG dominant". Figure |7^ shows the recovered F500 
values versus the model expectation for this subsample. The 
PZancA; measurements of the sample of 13, 104 MaxBCG clus- 
ters (not a subsample) is included for reference. 

Recently a new measure of cluster richness was developed 
with le ss scatter than the me asure presented in iKoester et alJ 
(I2007ah (iRykoffet all 1201 Ih . We test whether using this 
new richness measure will yield differing results with regard 
to measured versus expected SZ signal. A catalog with a 
new richness measure assigned to each MaxBCG cluster is 
available online|3 While for the previous richness measure 
the cluster richness-mass relation was calibrated using weak 
lensing, this new rel ation was calibrate d with an abundance 
matching technique (iRykoff et al.ll201 11) . Calibration of this 
new measure via weak lensing is still in progress. We use the 
richness-mass relation given in Eq. B6 of lRvkoffetai] dMT]) 
to determine Msooc and subsequently /?5ooc- Figure |7]3 shows 
the measured and expected Y^qq values using this new rich- 
ness measure for the 474 MaxBCG clusters. The red triangles 
show the Planck measurements using the old richness mea- 

clusters are only in Stripe 82, and thus are located towards the central part of 
the ACT map which has lower noise. 
http://kipac.stanford.edu/maxbcg 
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Fig. 7. — Left (a): Shown are the recovered Yjoo values as in Figure|5] except using a subsample of 126 "BCG dominant" clusters from the 474 MaxBCG 
clusters used above. Both the measured values (black circles) and expected values (blue squares) change using the new subsample as compared to Figure|5] The 
red triangles are the measured values from Planck for the 13,104 MaxBCG sample. Right (h): Shown are the same points as in Figure|5] however, using the new 
richness measure for the MaxBCG cluster sample given in Rvkoff et aljj2011i) . Note that both the measured values (black circles) and expected values (blue 
squares) change using the new richness measure as compared to Figure |5] The red triangles are the Planck measured values using the old richness measure for 
reference. 



sure for reference. Note that Planck would have different re- 
sults if they use the new richness measure. 

5.4. Contamination from Infrared and Radio Galaxies 

To investigate whether infrared galaxies may be reducing 
the SZ decrement at the MaxBCG cluster positions, we re- 
cover the F500 values from the 474 MaxBCG clusters studied 
above using the single band 218 GHz ACT map alone. We 
compare that to Y^qq values extracted at 474 random positions 
within the 218 GHz map PI A positive correlation between 
MaxBCG clusters and infrared galaxies would result in neg- 
ative 1500 values compared to the random sample. We find 
that both the MaxBCG cluster sample and the random sample 
have an SZ flux consistent with zero in the 218 GHz map, and 
we do not detect any significant excess of infrared signal cor- 
related with the MaxBCG sample. This is shown in Figure |8] 

We also cross-correlate the MaxBCG cluster catalog with 
the VLA FIRST catalog of radio sources to investigate how 
much radio galaxies may be reducing the measured SZ decre- 
mentH The FIRST survey uses the NRAO Very Large Array 
and covers over 10,000 square degrees of sky to a sensitivity 
of about 1 mJy at 1 .4 GHz. This survey also overlaps with the 
Sloan Digital Sky Survey. We cross-correlate to find the frac- 
tion of 474 MaxBCG clusters that have a radio source above 
a given flux threshold, within 5' of the identified cluster po- 
sition. Ass uming a typical spectral index for radio sources of 
-0.7 (e.g., Condonlll984t iLin et al.ll2009h . we choose a flux 
threshold cut of 50 mJy at 1.4 GHz to yield sources above 
2 mJy at 150 GHz. Such a source convolved with the ACT 
beam would have a temperature increment of about 3Q^K at 

^' For these measurements of Yfoo, we do not divide out the amplitude of 
the frequency dependence of the SZ signal, which is close to zero at the null 
frequency. So these are really measurements of — Ar/Tcmb. 

http://sundog.stsci.edu/first/catalogs/readme_12febl6.html 



148 GHz. Since the typical SZ signal from a cluster is about 
100/iA' to within a factor of a few when smoothed with the 
ACT beam, a 2 mJy radio source would start to significantly 
reduce the SZ decrement. We find that about 10% of the 474 
MaxBCG clusters investigated above have such a radio source 
within 5' of its identified center. This small correlation is not 
enough to explain the large discrepancy between measured 
and expected SZ signals shown in Figure |5] although it may 
be a contributing factor 

6. SIMULATED PLANCK SZ SIGNALS 
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Fig. 8. — Measured Fsoo values from 474 MaxBCG clusters using the 218 
GHz ACT equatorial map alone (black circles). For comparison is shown the 
same 7500 recovery procedure performed at 474 random locations in the 218 
GHz map (purple squares). 
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Fig. 9. — Recovered Y^qq values from 984 simulated clusters embedded 
in simulated CMB maps at 148 and 218 GHz with Planck-\ike noise (green 
diamonds) or the corresponding level of white-noise (purple hexagons). The 
maps were convolved with ACT beams as well as a 5' Gaussian beam, so that 
they correspond to Planck resolution. The recovered positions are offset from 
the true positions with a uniform random distribution between and 1.5 Mpc. 
The input yjoo values are shown as blue squares, and red triangles are the 
measure d values from the Planck satellite for a sample of 13,104 MaxBCG 
clusters JPlanck Collaboration et al.l2011d) . 



To investigate what Planck's measured SZ signal would 
be if there existed the amount of offset modeled in Fig- 
ure |6] we again use simulations. We embed 492 simulated 
MaxBCG clusters in two sets of simulated CMB maps at 148 
and 218 GHz that have Planck-like instrument noise added 
to themH We model the Planck noise using the the noise 
power spectra at 143 and 2 1 7 GH z shown in Figure 35 of 
iPlanck HFI Core Team etall (1201 Ih . We allow for an offset 
between the cluster SZ peaks and the identified cluster cen- 
ters that has a uniform random distribution between and 1 .5 
Mpc, analogous to Figure|6] We also convolve both maps with 
ACT beams and a 5' Gaussian beam to approximate Planck 
resolution. The results of the extracted isoo values are shown 
as green diamonds in Figure |9] We also show the case where 
the simulated clusters are embedded in simulated CMB plus 
white noise maps, using white noise levels that are similar 
to Planck noise levels. These results are shown as purple 
hexagons in Figure |9] 

Planck-liks noise is nearly white at these frequencies, so 
there is not much difference between the Planck-like case and 
the white noise case in Figure|9] The results of the Planck-like 
case in Figure|9]are also very different from those of the sim- 
ulated ACT case shown by the open black circles in Figure |6] 
The two differences between the simulations are the beam and 
the noise. In the ACT data there is 1/f noise, atmospheric 
noise, and noise from the primary microwave background^ 
All of this results in a redder noise spectrum than Planck's. 
The presence of red noise in the maps causes the matched 
filter to suppress more power on large scales than would be 
the case for white noise. This causes the filter to return a 
lower signal than actually exists if the signal is extracted from 

^' We doubled the number of maps and thus the cluster sample to shrink 
the error bars. This results in 984 embedded clusters in total. 

The instrumental noise in the ACT 218 GHz map is not low enough to 
remove all the primary microwave background signal from the 148 GHz map. 



a position that is offset from the cluster center If the signal 
is extracted at the cluster center, however, the matched filter 
will return the correct signal regardless of whether the noise 
is white or somewhat red (as can be inferred from Figure [TlFl 

We demonstrate this effect in Figure [10] where the blue 
solid curve is a 5' Gaussian profile. This is roughly the fil- 
ter response for a cluster profile convolved with the Planck 
beam as Planck's noise is close to white (see Eq.[T]i- The red 
dashed curve shows a 1 .4' Gaussian profile, which had been 
Fourier transformed and had all Fourier modes set to zero for 
I < 2000, roughly analogous to the effect of the ACT noise. 
The profile was then Fourier transformed back. One can see 
how a low flux recovery results from an offset from the center 
The green dotted curve shows the same as the red curve except 
for a 5' Gaussian profile. The low signal from a miscentered 
position is still apparent. This is why merely smoothing the 
ACT maps to match the Planck beam would not allow us to 
reproduce the Planck measurement. 

It is interesting to see from Figure |9] that the amount of 
offset modeled in Figure |6] to match the ACT data, can also 
explain the discrepancy found by Planck. If this amount of 
miscentering actually exists, then it would be the sole expla- 
nation of the discrepancies. However, this amount of offset 
between BCGs and SZ peaks is much larger than the distri- 
bution shown in Figure [2] So it may be that the SZ signal 
is intrinsically low by some amount or that the optical weak- 
lensing mass calibration is biased high. It is also possible that 
the fraction of false detections in the optically-selected sam- 
ple is larger than 10%, or some non-zero amount of radio or 
infrared galaxy contamination is filling in the SZ decrements 
(as discussed in Section l54l i. If so, these effects would serve 




Offset (arcmin) 

Fig. 10. — Shown is the filter response for differing beams and noise mod- 
els. The blue solid curve is a 5' Gaussian profile, which roughly approximates 
the filter response for a cluster profile convolved with the Planck beam since 
Planck's noise is close to white (see Eq.[T). The red dashed curve shows a 
1.4' Gaussian profile, which has been Fourier transformed and has had all 
Fourier modes set to zero for / < 2000, roughly analogous to the effect of the 
ACT noise. The profile was then Fourier transformed back. The green dotted 
curve shows the same as the red curve except for a 5 ' Gaussian profile. 

Note that estimating the noise power spectrum as CMB plus white noise 
when filtering the ACT data results in very large eiTorbars since the matched 
filter is no longer optimal. High-pass filtering the ACT maps enough to mini- 
mize the effect of the red noise causes ringing of the true cluster signal, which 
also results in a low measurement if the signal is exfi'acted from a position 
offset from the center. 



ACT Optical Richness - SZ Relation 



9 



to lower the measured SZ signal by the same amount for both 
ACT and Planck. In that case, a smaller amount of miscen- 
tering would be required to match the remaining discrepancy 
found by these instruments. It is thus likely that a combina- 
tion of effects may be at work resulting in the measurements 
found by ACT and Planck. 

In this vein, we have checked with simulations whether an 
intrinsic SZ signal as low as measured by Planck, combined 
with the distribution of BCG/SZ peak separations shown in 
Figure |2] can roughly match the SZ signal measured by ACT. 
We did this by lowering the SZ signal of the simulated clusters 
until they matched the Planck measurements, and then recov- 
ered this signal using wrong gas peak positions given by the 
distribution in Figure |2] as done in Section l4~2l We find that 
the recovered simulated SZ signal is at least a factor of two 
higher than measured by ACT. 

7. DISCUSSION 

From the analysis of simulated clusters embedded in ACT 
data presented in Figure [T] we expect robust recovery of the 
SZ flux from optically-selected MaxBCG clusters. This is as- 
suming that the position of the center of the gas concentra- 
tion is known. From the analysis of 52 MCXC clusters that 
fall within the ACT Equatorial and Southern survey regions 
shown in Figure ID we find agreement between the expected 
and measured Ysqqc-Mso qc relation. Such agreement is con - 
sistent with that found in Planck Collaboration et al.l ( 1201 Ib l). 
Both of these figures taken together give confidence in the SZ 
flux recovery pipeline and the analysis of the ACT maps. 

In Figure |5] we find that the recovered ACT SZ flux from 
474 optically-selected MaxBCG clusters is lower than both 
th e model e xpectatio ns and the measu red Planck values given 
in lPlanck"C ollaborat ion et al.l (1201 Id) . Since we expect from 
Figure |2] some offset between the positions of the BCGs and 
the centers of the gas concentrations, we explore this possibil- 
ity in more detail. Given the difference in resolution and noise 
properties between ACT and Planck, such an offset would re- 
sult in different measured SZ signals between the two instru- 
ments (see Figure [Toll. 

Figure|2]gives the offset distribution for clusters in common 
between the fufl MaxBCG and MCXC cluster catalogs. Using 
simulations. Figure |3] shows that this amount of offset is not 
enough to explain the measu red SZ signal shown in Figure |5] 
Modeling the offset using the' Johnston et al.l (12007 b) distribu- 
tion, which only includes misidentification of the BCG by the 
optical cluster finding algorithm, produces a slightly smaller 
amount of discrepancy between expected and measured SZ 
signals than what is shown in Figure |3] To match the mea- 
sured SZ signal found by ACT requires an offset distribution 
that is significantly larger A distribution that is uniformly ran- 
dom between and 1 .5 Mpc gives a better fit to the ACT data, 
as shown in Figure |6] This amount of offset also can explain 
the Planck measured discrepancy as shown in Figure |9] 

It is possible that the subsample of clusters in common be- 
tween MCXC and MaxBCG catalogs is not representative 
of the full MaxBCG cluster sample. This subsample is in 
fact special in that these clusters are found using both optical 
and X-ray selection techniques. When Planck measured the 
SZ signal for such a subsample, no significant discrepancy 
between the measured and expec ted SZ signals was found 
dPlanck Collaboration et al.l201 Id) . This is in contrast to their 
results for the full optically-selected sample. This suggests 
that the subsample of clusters in common to both catalogs is 
not representative of the full optically-selected sample, at least 



in some regard. 

We also assume throughout this work that the lAmaud et al.l 
( 120101) profile accurately describes the gas profiles of clusters 
in the MCXC and MaxBCG catalo gs. The good agreement we 
find in Figured and also shown in Planc k Coflaboration et al.l 
(1201 lb) suggests that this profile works well for clusters in 
the MCXC sample. However, this profile may not hold for 
the optically-selected MaxBCG sample. In particular, to mea- 
sure a lower SZ signal with ACT than found by Planck would 
require that the gas distribution is wider (more spread out) 
than that given by the Arnaud et al. (2010) profile. In such a 
case, given the finer resolution of ACT, it is conceivable that 
ACT could measure a lower signal than Planck when the SZ 
signal is extracted using an incorrect narrower profile for b oth 
datasets. However, since we expect clusters to be in hydro- 
static equilibrium to first order, the consequence of a wider 
gas distribution is a lower total mass within Msoot- if the gas 
mass fraction is fixed. This would imply lower normalizations 
of the Lx-MsQOc and A^200m - A^sooe relations. That in turn 
would result in a lower gg value obtained from the MaxBCG 
sample than reported in lRozo et al.l (12010 ). the latter of which 
is currently in agreement with erg values obtained from X-ray- 
selecte d cluster sa mples (Henry et al. 2009; Vikhlinin et al] 
2009: Mantz et alJiSoU) . 

Some other recent papers have added more data and anal- 
ysis shedding ligh t on the measured d i screpa ncy in SZ signal 
found by Planck. iBiesiadzinski et all ( 120 12h argue that there 
is not really a discrepancy if on e considers t he 2a uncertainty 
on the N2oom —Msooc relation in lRozo et al.l ((2009). However, 
if Rozo et al.. ( 2010.) used an A^200m -^500c relation that had a 
normalization that is 2(7 away from the true value, then again 
it is surprising that they find a constraint on with an op- 
tical cluster sample that is in such good agreement with the 
(78 constraints found using X-ray-selected samples. Such a 
shift in the A^200m - A^sooc normalization would also put the 
derived cos mological constrain ts in some tension with previ- 
ous results (lRozoetal.ll2012ah . lAngulo et aTl ( 120121) . using 
Millennium-XXL s i mulat ions, argue along the same lines as 
IBiesiadzinski et all (120121) suggesting that the normalization 
of the richness-mass relation is incorrect (see last paragraph 
of Section 4 in that w ork). 

Bauer et aH (120121) tested the A^200/>i —Msooc relation by veri- 
fying the weak-lensing cluster mass estimates upon which this 
relation was calibrated. They did this by using gravitational 
lensing magnification of type I quasars in the backgrounds of 
these clusters. T heir results support the mass normalization 
of thelRozo et al. , (20 09) N200m—Msmc relation. Recent work 
by lRozo et alJ (l2012d) suggests that part of the Planck mea- 
sured discrepancy may be due to different X-ray instrument 
calibrations between Chandra and XMM-Newton and/or dif- 
fering analysis procedures when dealing with data from each 
instrument (see Figure 3 of that work). Rozo et al. (2012d!) 
argue that X-ray scaling relations derived from Chandra ob- 
servations may be systematically high compared to those from 
XMM-Newton observations. Since the A^zoom - sooc relation 
was calibrated using an Lx-Msfxk- prior derive d from Chan- 
dra observations (Vikhlinin et al. 2009), and the lArnaud et al.l 
(2010) profile was derived using XMM-Newton observations, 
this may be responsible for part of the discrepancy. Rozo et al] 
(2012a,b,c) include further discussion in this direction and 
claim that three ingredients may be responsible for the ten- 
sion between measured and expected SZ signals measured 
by Planck for the MaxBCG sample. These ingredients are a 
mass-richness relation that is normalized high by 10%, Chan- 
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dm derived X-ray masses from lVikhlinin et al.l (12009b that are 
biased low by 20% due to non-thermal pressure, and XMM- 
Newton X-ray derived masses that are ~ 10 to 20% lower than 
those derived from Chandra. More data will clarify this issue. 

Given the above discussion, the interpretation of both the 
ACT results presented here and the Planck results stacking 
MaxBCG clusters remains unclear Since each possibility 
mentioned above, on its own, would be unlikely to explain 
the full discrepancy between expectation and measurement of 
both data sets, it is likely that multiple effects are responsi- 
ble for the low SZ signals. It is possible that a combination 
of factors is lowering the SZ signal equally for both ACT and 
Planck, and that offsets between gas peaks and BCGs result in 
the remaining discrepancy between ACT and Planck measure- 
ments. We have shown that even in this scenario however, a 
larger offset distribution than is seen in X-ray-selected cluster 
samples is still required to explain the discrepancy between 
ACT and Planck measured signals. 

Further investigation is required using multi-wavelength 
techniques to fully contextualize the SZ signal measured by 
both Planck and ACT for this sample of optically-selected 
clusters. High-resolution SZ instruments such as CARMA 
and MUSTANG2 as well as more sensitive datasets such as 
those forthcoming from ACTPol and SPTPol will provide 
more clarity. This will yield a better understanding of the 
astrophysics of galaxy clusters. Such knowledge will also be 
vital to assess the use of galaxy clusters as probes of structure 
growth and the cosmological parameters of the Universe. 

NS would like to thank James Bartlett, Sarah Hansen, Jean- 



Baptiste Melin, Eduardo Rozo, Kendrick Smith and Risa 
Wechsler for useful discussions as well as Rocco Piffarett and 
Eli Rykoff for providing the MCXC catalog and the new rich- 
ness MaxBCG catalog respectively. NS would also like to 
thank the organizers of the KITP Monsters Inc. conference, 
which facilitated further discussions. 

This work was supported by the U.S. National Sci- 
ence Foundation through awai'ds AST-0408698 for the ACT 
project, and PHY-0355328, AST-0707731 and PlRE-0507768 
(award number OlSE-0530095). The PIRE program made 
possible exchanges between Chile, South Africa, Spain and 
the U.S. that enabled this research program. Funding was also 
provided by Princeton University, the University of Pennsyl- 
vania, and a Canada Foundation for Innovation (CFl) award to 
UBC. Computations were performed on the GPC supercom- 
puter at the SciNet HPC Consortium. SciNet is funded by: the 
Canada Foundation for Innovation under the auspices of Com- 
pute Canada; the Government of Ontario; Ontario Research 
Fund - Research Excellence; and the University of Toronto. 
ACT is on the Chajnantor Science preserve, which was made 
possible by the Chilean Comision Nacional de Investigacion 
Cientifica y Tecnologica (CONICYT). 

NS is supported by the National Science Foundation under 
Award No. 1 102762. During the completion of this work, NS 
was also supported by the U.S. Department of Energy contract 
to SLAC no. DE-AC3-76SF005 15 and in part by the National 
Science Foundation under Grant No. 1066293 and the hospi- 
tality of the Aspen Center for Physics. The data will be made 
public through LAMBDA (http://lambda.gsfc.nasa.gov/) and 
the ACT website (http://www.physics.princeton.edu/act/). 



REFERENCES 



Andersson, K. et al. 2011, ApJ, 738, 48, 1006.3068 

Angulo, R. E., Springel. V., White, S. D. M., Jenkins, A., Baugh, C. M., & 

Frenk, C. S. 2012, ArXiv e-prints, 1203.3216 
Amaud, M., Pratt, G. W., Piffaretti, R., Bohringer, H., Croston, J. H., & 

Pointecouteau, E. 2010, A&A, 517, A92, 0910.1234 
Bauer, A. H., Baltay, C, Ellman, N., Jerke, J., Rabinowitz, D., & Scalzo, R. 

2012, ArXiv e-piints, 1202.1371 
Benson, B. A. et al. 2011, ArXiv e-prints, 1112.5435 

Biesiadzinski, T., McMahon, J. J., Miller, C. J., Nord, B., & Shaw, L. 2012, 

ArXiv e-prints, 1201.1282 
Bohringer, H. et al. 2007, A&A, 469, 363, arXiv:astro-p h/0703553| 
Bonamente, M. et al. 2011, ArXiv e-piints, 1112.1399 
Bonamente, M., Joy, M., LaRoque, S. J., Carlstrom, J. E., Nagai, D., & 

Marrone, D. P 2008, ApJ, 675, 106, 0708.0815 
Carlstrom. J. E., Holder, G. R, & Reese. E. D. 2002, ARA&A, 40, 643, 

arXiv:astro-ph/0208192 
Condon, J. J. 19S4, ApJ, 287, 461 
Das, S. et al. 2011, ApJ, 729, 62, 1009.0847 

Draper, P, Dodelson, S., Hao, J., & Rozo, E. 201 1, ArXiv e-prints, 

1106.2185 
Dunner, R. et al. 2012, in prep. 
Fowler, J. W. et al. 2010, ApJ, 722, 1148, 1001.2934 

. 2007, Appl. Opt., 46, 3444, arXiv:astro-ph/0701020 

Haehnelt, M. G., & Tegmark, M. lyyb, MJNKAS, 2/9, 543, 

arXiv:astro-ph/9507077 
Haiman, Z,., Mohr, J. J., & Holder, G. P 2001, ApJ, 553, 545, 

arXiv:astro-ph/0002336 
Hajian, A. et al. 2U11, ApJ, 740, 86, 1009.0777 
Hall, N. R. et al. 2010, ApJ, 718, 632, 0912.4315 
Hand, N. et al. 2011, ApJ, 736, 39, 1101.1951 

Henry, J. P, Evrard, A. E., Hoekstra, H., Babul, A., & Mahdavi, A. 2009, 

ApJ, 691, 1307, 0809.3832 
Holder. G., Haiman, Z., & Mohr, J. J. 2001, ApJ, 560, Llll, 

arXiv:astro-ph/0105396 
Johnston, L). E., Sheldon, E. S., Tasitsiomi, A., Frieman, J. A., Wechsler, 

R. H., & McKay, T. A. 2007a, ApJ, 656, 27, arXiv:astr o-ph/0507467, 
Johnston, D. E. et al. 2007b, ArXiv e-prints, 07uy. 1 139 " 
Koester, B. P et al. 2007a, ApJ, 660, 239, arXiv:astro-ph/0701265| 

. 2007b, ApJ, 660, 221, arXiv:astro-pn/0/012ti8 

Komatsu, E. et al. 2011, Apjs, 192, IS, 1U01.43JS 

Lin, Y, Partridge, B., Pober, J. C, Bouchefry, K. E., Burke, S., Klein, J. N., 
Coish, J. W., & Huffenberger, K. M. 2009, ApJ, 694, 992, 0805.1750 



Lueker, M. et al. 2010, ApJ, 719, 1045, 0912.4317 

Mandelbaum, R., Seljak, U., & Hirata, C. M. 2008a, J. Cosmology 

Astropart. Phys., 8, 6, 0805.2552 
Mandelbaum, R. et al. 2008b, MNRAS, 386, 781, 0709.1692 
Mantz, A., Allen, S. W., Rapetti, D., & Ebehng, H. 2010, MNRAS, 1029, 

0909.3098 

Marriage, T. A. et al. 2011, ApJ, 737, 61, 1010.1065 

Marrone, D. P et al. 201 1, AiXiv e-prints, 1 107.51 15 

Melin, J.-B., Bartlett, J. G., & Delabrouille, J. 2006, A&A, 459, 341, 

arXiv:astro-ph/0602424 
Melm, J.-Ji., Bartlett, J. U., Delabrouille, J., Amaud, M., Piffaretti, R., & 

Pratt, G. W. 2011, A&A, 525, A139, 1001.0871 
Motl, P M., Hallman, E. J., Bums, J. O., & Norman, M. L. 2005, ApJ, 623, 

L63, arXiv:astro-ph/0502226 
Nagai, U. 20Ub, ApJ, b3U, 3jS, arXiv:astro-ph/05 12208 
Piffaretti, R., Amaud, M., Pratt, U. W., fointecouteau, b., & Melin, J.-B. 

2011, A&A, 534, A109, 1007.1916 
Planck Collaboration et al. 201 la, A&A, 536, Al 1, 1 101.2026 

. 2011b, A&A, 536, AlO, 1101.2043 

. 2011c, A&A, 536, A12, 1101.2027 

Planck HFI Core Team et al. 201 1, A&A, 536, A6, 1 101.2048 

Pratt, G. W., Croston, J. H., Arnaud, M., & Bohringer, H. 2009, A&A, 498, 

361, 0809.3784 

Rozo, E., Bartlett, J. G., Evrard, A. E., & Rykoff, E. S. 2012a, ArXiv 

e-prints, 1204.6305 
Rozo, E., Evrard, A. E., Rykoff, E. S., & Bartlett, J. G. 2012b, Ai'Xiv 

e-prints, 1204.6292 
Rozo, E., Rykoff, E. S., Bartlett, J. G., & Evrard, A. E. 2012c, Ai'Xiv 

e-prints, 1204.6301 
Rozo, E. et al. 2009, ApJ, 699, 768, 0809.2794 

Rozo, E., Vikhlinin, A., & More, S. 2012d, ArXiv e-prints, 1202.2150 

Rozo, E. et al. 2010, ApJ, 708, 645, 0902.3702 

Rykoff, E. S. et al. 2011, ArXiv e-prints, 1104.2089 

Sehgal, N. et al. 2011, ApJ, 732, 44, 1010.1025 

Sheldon, E. S. et al. 2009, ApJ, 703, 2217, 0709.1153 

Sifon, C. et al. 2012, ArXiv e-prints, 1201.0991 

Swetz, D. S. et al. 2011, ApJS, 194, 41, 1007.0290 

Vikhlinin, A. et al. 2009, ApJ. 692, 1060, 0812.2720 . 

Voges, W. et al. 1999, A&A, 349, 389, arXiv:astro-ph/99()93T5l 

Wang, L., & Steinhardt, P J. 1998, ApJ, 30S, 4SJ, arAiv:astro-ph/9804015 1 
York, D. G. et al. 2000, AJ, 120, 1579, arXiv:astro-ph/u006396i 



